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FOREWORD 

This report summarizes studies performed under Contract NAS12-500 with 
the Guidance Laboratory, NASA Electronics Research Center, Cambridge, 
Massachusetts. The period of performance was from 15 September 1966 through 
31 October 1967. Technical director for the contract was Mr. W. E. Miner of 
the Electronics Research Center Guidance Laboratory. 

Nortronics personnel who contributed to this study and this report are 
listed on the cover page. The writing and organization of this report was 
done by M. L. Thompson, 

Computer time needed during the course of this study was obtained from 
the Marshall Space Flight Center's Computation Laboratory under the NASA 
Resource Sharing Program. 
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SUMMARY 

Described in this report are studies completed under Contract NAS12-500 
with the NASA Electronics Research Center, Cambridge, Massachusetts. 

These studies are concerned with series solutions of two-point boundary 
condition problems that result from applying the calculus of variations to 
optimal guidance problems. The major objective was to use these solutions to 
devise accurate, analytical approximations to guidance functions for optimal 
ascent to orbit. This work is essentially an extension of that completed 
under Contract NASW-1165, which is described in reference 1. 

An outstanding problem associated with these studies is the need to perform 
algebraic manipulations of lengthy and complicated formulas . In particular 
the multiplication of symbolic, polynomial- like expressions occurs frequently. 

An integral part of the work was to devise new techniques or apply existing 
ones to automate symbolic mathematical operations on digital computers. 

In addition, research in the calculus of variations was carried out by 
Mr. Robert Silber. The result of these efforts was a modification of the 
Bolza problem to include end-orbits governed by n-body attractions. This work 
is available in reference 2 and is not discussed at length in this report. 
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SECTION I 

GENERAL STATEMENT OF GUIDANCE PROBLEM AND APPROACH TO SOLUTION 

In this section, only a cursory description of the problem treated and 
the approach to solution is given. A detailed description has been given in 
reference 1, 

1.1 DEFINITIONS 

We make the following assumptions and definitions regarding the vehicle 
and its performance: 

a Constant thrust magnitude, F, and mass flow rate, m. 

• Rocket is represented by a point mass, m. 

• Single-stage rocket with continuous thrusting to orbit. 

• Spherical, non- rotating, airless central body of attraction. 

• Trajectory and orbit are coplanar. 

• The optimum trajectory is the one that requires a minimum time from 
initial state to- terminal orbit. 

1.2 DIFFERENTIAL EQUATIONS 

The coordinate system used is shown in the following figure. The equations 
of motion are written in this system 
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•• F . 

X - sxn X - 


y = - cos X " g,, 
m y 


Application of the calculus of variations yields an additional system of 
differential equations. These equations may be coupled with the preceding 
equations to yield; in vector form: 


7.2 ~ ■- 

” m 1 


A = + (A • R)y R 


( 1 ) 


where 



The second of equations (1) is usually called the Euler-Lagrange equation, and 
the A and A are called Lagrange multipliers . 


1.3 BOUNDARY CONDITIONS 

The initial conditions on the solutions of equations (1) are known with 
the exception of l(t^) and A(t^). At the unspecified terminal time, t^, it is 
necessary that the solution of equations (1) satisfy functions that generally are 
of the form 

F^(R, V, A,' I) =0 (j - 1,...,4) (2) 

Some of these are the transversality conditions and some are relations 
that describe the desired terminal orbit. 
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1.4 APPROACH TO SOLUTION 

We first note that there are as many boundary conditions to be satisfied 
as there are unknown initial conditions . Because of the homogeneity of the 
Euler-Lagrange equation, one of the multipliers may be set to an arbitrary 
value at the initial point. The same thing can be accomplished by letting 
jxj = 1 at t = t^. This additional condition will also introduce simplifications 
in subsequent developments. For brevity we will refer to this relation as the 
"scaling condition". 


We now have a total of five conditions on the solutions of equations (1) 

and also five unknowns. In order to make use of the we expand them in 

Taylor series about the interval (t. - t ), to express them approximately at 

f o 

the initial point. This yields a system of equations of the form: 


k 


I 

n = 0 


i- 

n! 



(F.) 



^0^ 


n _ 


(3) 


Repeated substitution of equations (1) and their derivatives in the series 

coefficients of equations (3) is used to reduce equations (3) to a system of 

algebraic equations in the Lagrange multipliers. Equations (3), together with 

the scaling condition, then implicitly define values for the multipliers and 

the flight time, (t^ - t ), which we define as At. 

t o 


Numerical verification of the validity of this approach for ascent to 
circular orbits has been demonstrated previously for numerous trajectories 
(see reference 1) . Usually the error in the solution for the multipliers is 
1 percent or less for series of order four or greater. 
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Actually , the solution for the multipliers was not carried out simultan- 
eously with the solution for At. Instead, one of the series expansions was 
solved for At explicitly by a procedure known as "series inversion" and the 
result was substituted into the other expansions. The reason for this is dis- 
cussed in the next section. 

Now, with At eliminated from the equations we have four algebraic equations 
in the four unknown Lagrange multipliers. These equations can be solved for 
values of the multipliers at the initial point when values of the state 
variables and constants are known. It can also be shown that a simple relation- 
ship exists between the multipliers and the steering angle, x fot an optimal 
trajectory [cf. equations (1)]. Thus,- a solution for the multipliers is a 
solution to the optimal guidance problem. 

Subsequent sections of this report discuss the applicatipn of this approach 
to three different missions; optimal ascent to circular orbit, optimal ascent 
to rendezvous with a circular orbiting. target satellite, and optimal ascent 
to elliptical orbit. 

In order to test or verify the validity of the analytical developments in 
this study, comparison is made with numerically computed, optimum trajectories. 
These trajectories are called "nominal". Values for the state variables and 
constants are taken from these nominal trajectories and used in the numerical 
evaluation of formulas . Comparisons are then made between nominal values for 
the multipliers and flight time and those values predicted by the formulas. 
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SECTION II 

OPTIMAL ASCENT TO CIRCULAR ORBIT 


2 . 1 INTRODUCTION 

In this guidance problem we are concerned with the problem of steering 
a vehicle along an optimum path to insertion into a circular orbit of specified 
altitude. 


The terminal functions to be satisfied are; 


,,2 2 

- V 

f CO 


= 0 


2 2 

R^ - R 

f CO 


= .0 


(4) 


R, 



= 0 


+ A^y ~ = 0 


It can be shown that the last of equations (4), the transversality condition, 

is for this problem, an Integral of the differential equations and is valid at 

to* It can be applied directly without being expanded in a series, as can the 

scaling condition, |a| = 1, by definition at t = t . 

o 


The first of equations (4) will be used to illustrate the expansion in a 
Taylor series ■ 
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= 0 = (V^ - ) + 2(J-Y) (At) 

f CO CO O ' o 


+ (Y-V + V*V) (At)^ + ~ (Y-V + 3V*V) ) (At)"^ + ... (5) 

o j ! o 


where At = (t^ - t ) and the subscript "o" indicates a value at t = t . 

After the remaining functions are expanded in a similar fashion, the 
derivatives of the differential equations can be substituted into the series'! 
coefficients Co reduce equations (4) to a system of five nonlinear algebraic 
equations in the four Lagrange multipliers and At. (For reference, these 
derivatives are listed in Appendix A.) 

At this point we could solve these five simultaneous equations. However, 
we choose to eliminate At from. the equations by a technique called "series 
inversion". One of the series is inverted to obtain an analytical expression 
for At and this expression is then substituted for At in the other expansions. 
There are two reasons for this approach. First, one of the series may, for a 
given order, implicitly define more accurate values for At than the other 
series. In fact thds was show to be the case in reference 1. Second, it may 
be desired to calculate approximate At values independently of the multipliers 
during the guidance process. 


2.2 SOLUTION FOR MULTIPLIERS 

The solution of the systems of algebraic equations for the four Lagrange 
multipliers was discussed at length in reference 1. Details to be found there 
will not be discussed here. 


2-2 



TR-792-8-283 


NORTRONICS - HUNTSVILLE 


Of major concern during this study was the possible simplification of 
analytical derivations. It has been verified previously that systems of 
equations can be derived that define quite accurate solutions for the multipliers 
More accuracy can be attained by using higher order series. Therefore, primary 
attention was directed toward analytical solutions of the equations that define 
the multipliers, elimination of insignificant terms, and determining approxi- 
mations associated with certain classes of missions . For convenience and 
simplicity, the systems of equations derived from third-order series were used 
in the search for analytical solutions. • 

Two methods for solution of the equations for the multipliers were examined 
The method of "Successive Substitutions" was considered as a means of getting 
an explicit formula for the multipliers as functions of the local state 
variables and mission constants. At first it seemed that this technique could 
be easily applied. But in order to obtain the third step approximations, about 
100,000 terms would be involved. A few algebraic simplifications such as 
combining terms were possible but did not cause any significant reduction in 
numbers of terms. It was decided that, if used, this method should be applied 
after the coefficients in the equations were numerically evaluated. 

The second method considered was to invert the system of equations in the 
multipliers. This would lead to less complicated formulas for the solution 
but not much is known about the practical question of convergence. Formulas 
for this inversion are derived in Appendix B. For illustrative purposes, the 
formulas are for two series in two variables. The inversion technique should 
converge rapidly if estimates of the multipliers are available. By restricting 
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the problem to a class of trajectories it may be possible to expand the 
equations in the multipliers about a "standard estimate" associated with a 
class of trajectories. One would then have a system of equations in which the 
variables are corrections to the estimates. If the corrections are relatively 
smaller than the estimates, there should be many opportunities for simplification 
and approximation. 

An approach analogous to the above (but in one variable) was taken in the 
solution for time-to-go or At, described in the following subsection. There, a 
great amount of simplification was possible. 

2.3 SOLUTION FOR FLIGHT TIME (At) 

Of the three possible expansions to be inverted for At, the one describing 
the terminal velocity condition, equation (5) was chosen. The reasons for this 
choice were: (1) A third-order series implicitly defines accurate At values. 

(2) Velocity, for an optimal trajectory, should be a monotone increasing 
function of time and the possibility of a singularity is avoided. 

We rewrite equation (5) as 
Z = A^At + 

where 

Z = [v^ - V^I/2? • V 

\ CO ol 



A = (V • V + V ♦ V)/2V • V 

^ • ••• • • « • 

= (V • V + 3V * V)/6V * V. 
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Equation (6) can be inverted- to obtain an approximation for At 

At r B^Z + + B^Z^ + . . , (7) 

It often happens that the inverse series, equation (7), requires an 
unwieldy number of terms to approximate a root of equation (6); and the coef- 
ficients of the inverse series .become progressively more complicated. The 
number of terms needed in the inverse series can be significantly reduced by 
expanding equation (6), not about the origin but about an estimate of t^. 
Define 

At = T + dt, 

where x is an estimate and dt a correction to x. -Substitution into equation 
(6) yields 

Z = A^dt + dt^ + A^ dt^ (8) 

Equation (8) can now be solved for the correction by a series of inversion. 


There are several means of estimating At. One way is to take only the 
linear terras of equation (6) and have 

At = Z = T (A^ = 1). 

Another way is to use the -"rocket equation" which is fairly accurate for short- 
arc flights: 

m A. 

T - J <1 - B ), 

where 


A = -(''co - 

A third method of estimating time is to assume that it is known from a 
previous in-flight computation, say, five seconds ago, or that the expected 
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flight time is known from preflight planning. This latter possibility is 
quite attractive, because the At estimation formula is rather tolerant of 
"poor" estimates. 


Considerable simplification is possible in equation (8). It is only 
necessary, for most estimates, to invert the series to one term; i.e., 

t I 

Z = dt 


dt 


t 



■> 


where 

Z = Z - |a^t + ^2^^ ^3^^) 

' 2 
A^ = A^ + 2tA2 + 3t A^ 

' . . 2 

Within the A coefficient it was possible to neglect the term 3 t A and 

3 

frequently possible to neglect the term 2xk^. However, to assure accuracy 

over a reasonable range of T values , the latter term was retained with a further 

approximation. The formula for was approximated as 

k^~~V' V/2V * V = c 

The expression for At in simplified form is 

Z - (t + A + A T^l 

At = T + * — 

1 + 2C 

or 

At = T + dt., (9) 

A’ 
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In cases where It might be necessary to invert equation (8) to higher 

T 

orders than the first, the coefficients can be replaced with the coeffi- 
cients of equation (6) with no significant loss of accuracy. 

« 

The terms Z and C both involve terms with V in their denominators. This 
leads to formulas that involve A^^and the denominators when the equations 

of motion are substituted. To circumvent this difficulty, another approximation 
was made. It is assumed that the components of V, in denominators only, 
would be treated as measured accelerations rather than the accelerations 
associated with the optimal trajectory. This approximation was tested in 
equation (9) by using accelerations obtained by perturbations on the steering 
angle around a nominal, optimum value. The effect on At values was relatively 
insignificant as Table 2-1 indicates. 


Table 2-1. EFFECT ON At VALUES 


% ERROR IN X 

ERROR IN At (SEC) 

1 

0.45 

2 

0.85 

3 

1.5 

4 

1.6 

5 

1.66 

-1 

0.04' 

-2 

0.5 

-3 

1.1 

-4 

2.0 

-5 

2.7 
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No further simnl i.f icatf ons in equation (9) appear to be possible. For 
missions different from those represented by the nominal data used, the approxi- 
mations may be invalid. 

If equation (9) is written in terms of the multipliers after substitution 
of the differential equations, there is a total of 19 terms. (These are tabu- 
lated in Appendix C.) This is an improvement over results previously obtained 
(reference 1) where the At formula had 50 terms when written in terms of the 
multipliers . 

The great difficulty in eliminating At from the equations is substitution 
of a formula such as equation (9) into the other series where it must be 
squared, cubed, etc. But note that if t is a good estimate for At, then dt 

A 

is relatively smaller than t. We would then expect the following approximations 
to the powers of At to be sufficiently accurate. 

(At)^t 

(At)^ ~ +3x^dt. 

A 

etc. 

In fact, for a comfortable range of t, these formulas approximate the 
powers of At rather well. Representative data is shown in Figures 2-1 through 
2-3. 

It should be noted that At and its powers can be approximated by expressions 
that involve only 19 terms in the multipliers. Further, the algebraic form of 
the approximations is the same for the different powers . Only the coefficients 
need be slightly altered by a constant multiplicative factor. 
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The effect of these approximations on the solution for the multipliers 
was estimated without actually substituting the At formula. This was done by 
evaluating the approximations for .At and its powers with nominal data, including 
the multipliers. The coefficients of the equations in the multipliers were 
evaluated with these At values Instead of nominal At values, and then solved. 

The effect of the approximations for At was a small decrease in accuracy in the 
solution for the multipliers from that obtained when nominal At values were used. 
The difference between the two solutions ranged from less then one percent to 
about five percent in the worst case. 

2.4 SUBSTITUTION OF At EXPRESSION INTO OTHER SERIES 

As discussed previously, the expression for At is to be substituted into 
the other two series in order to have a system of four equations in the multi- 
pliers. From equation (3), it is evident that the coefficients of At in the 
series will generally be pol 3 momlal-like expressions in the multipliers. We 
are thus faced with the rather formidable problem of substituting one multi- 
variate polynomial into another. The problem can be represented as determining 
the right-hand side of the equation: 

Po + PfQi + P2Q2 + P3Q3 

where the P's and Q's are multivariable pol}nnomials . The right-hand side, S, 
should be simplified in the sense that like powers of the variables and like 
coefficients are collected. 

Note that the problem is considerably simplified if the approximations 
to At, previously described, are used. In that event, each of the Q's is 
essentially the same as far as algebraic manipulation is concerned. Only the 
coefficients associated with each Q need be distinguished . The problem is still 
not practically manageable if done by hand. 
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The Northrop MULPO program, developed under Contract NAS12-500, was used 
to perform these algebraic manipulations. (See Section V for a more detailed 
description.) The analytical approximations for At and its powers were sub- 
stituted into third-order expansions of the terminal functions that specify 
cutoff radius and path angle. The results are shown in Tables 2-2 and 2-3. (See 
subsection 5.3 for a key to reading Tables 2-2 and 2-3.) There is a total of 100 
different combinations of lambda exponents in one expansion (radius) and 176 
in the other (path angle) . As an estimate of the reduction in numbers of terms 
generated, the At expression without truncation of powers was substituted and 
the result simplified by MULPO. The resulting series had over 9000 terms. 

2.5 SIMPLIFICATION OF COEFFICIENTS 

Attempts were made to eliminate terras in the system of equations that 
define the initial values of the multipliers. (The formula for At had not been 
substituted.) This problem was discussed previously in reference 1. Essentially 
no simplifications other than those described in reference 1 were achieved. 

One approach that was tried was to write the differential equations' 
derivatives in vector form and then attempt to drop numerically insignificant 
terms. (See Appendix A.) This approach was less fruitful than working with 
components of the vectors as was' done previously. 

It was noted that simplifications that are made according to the be- 
havior on one trajectory are also valid for neighboring trajectories that 
have different initial values. Thus, simplifications already found are 
valid over a class or field of trajectories. 
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After the analytical approximations for At and. its powers are substituted, 
as shown in Tables 2-2 and 2-3, there may be possibilities for further simpli- 
fications. The system of equations in the multipliers obtained after substitution 
of At has not been coded for solution because of lack of available time. 

However, Northrop is sponsoring an effort to code a computer program to solve 
these equations numerically in order' to determine their validity. A part of 
this effort would also involve the elimination of terns that do not contribute 
significantly to the solution. 
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Table 2-2. SERIES FOR RADIUS CONDITION 
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... .JLff 

Off 

O.ff 

Off 

Off 

.Iff 
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Table 2-2. 

SERIES 

FOR RADIUS 

CONDITION (continued) 










LI 

L2 

L3 

L4 

L5 

L6 









C 

1 • 

0 

1 




1 













0 


0« 

0« 

0* 

Ilf 

1 ® 

0» 

1® 

0« 

Ok 

IK 






. . 





Ll 

L2 

’L3 

"L4 

L5 

L6 









0 

1 

I 

0 




1 













0 



0« 

1* . 

0» 


0® 

1® 

0« 

Ok 

IK 











Ll 

L2 

L3 

LA 

15 

L6 









0 

2 

0 

0 




1 













c 


0* 


o-» 

0 ^ 

1 ® 

0* 

1* 

0* 

Ok 

IK 




1 













D 


0* 

•1« ■ 

0* 

0* 

U 

o» 

Ilf 

6 * 

Ok 

iK 






’• * 


•' - 



U1 

L2 

L3 

LA 

15 

L6 









1 

0 

0 

1 




1 













D 


0* 


0« 

1* 


1* 

Oif 

0* 

Ok 

1 K 











Ll 

L2 

L3 

L'l 

45, 

It 









1 

0 

1 

0 




1 













D 


0» 


L* 


1* 

1 If 

0® 

Ok 

Ok 

iK 











Ll 

L2 

L3 

La 

Li; 

Lt. 









1 

1 

0 

0 




I 













C 


0* 

1« 

0» 

b» 


1« 

Ok 

Ok 

On 

IK 




1 













c 



0« 

o« 

0*:- 

1 » 

0« 

Ilf 

Ok 

Ok 

IK 




1 













D 


0» 

1« 

0« 

0* 

1® 

1® 

OK 

Ok 

Ok 

1 « 




1 













D 



0« 

0» 

0® 

i» 

0» 

iK 

Ok 

Ok 

IK 
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Table 2-2, SERIES FOR RADIUS CONDITION (continued) 










LI 

L2 

L3 

L4 

i. 5 

L6 









2 

0 

0 

0 




1 













c 

1 

1* ■ 


Os- 

b-r- 

! ft ' 

1ft 

Oft 

Oft 

Oft 

Ift 



D 


1» 

0« 

o» 

0 ft 

1ft 

1 ft 

Oft 

Oft 

Oft 

1ft 











LI- 

L2 

L3 

l4 

L5 

LA 









0 

• 1 

1 

I 




1 













B 

1 

0* 

0« 

0* 

Oft 

] ft 

o» 

1 ft 

1 ft 

1ft 

1ft 



C 

1 

0» 

0« 

0-R 

Oft 

1 ft 

Oft 

1ft 

1ft 

1ft 

1ft 



C 


0« 

0« 

0< 

Oft 

1 ft 

Oft 

1 ft 

1» 

1ft 

1ft 











LI 

L2 

L3 

LA 


le 









0 

1 

2 

0 




1 













B 

1 

0* 

0« 

■ 0«* 

Oft 

1 ft 

Oft 

1ft 

2ft 

Oft 

1ft 



C 

1 


0« 

0* 

Oft 

1ft 

0» 

1ft 

■ 2ft 

Oft 

l« 



D 


0» 

O^S 

Oft 

Oft 

1 ft 

Oft 

l^J 

2ft 

Oft 

1ft 






-- 





L'l 

L2 

L3 

1 A 

L? 

Lt 









0 

. 2 

0 

.. 1 




1 













0 


o«- 

2 « 

Oft 

1ft 

1 ft 

Oft 

0 ft 

0« 

Oft 

1ft 











LI 

L2 

L3 

LA 

L'S 

L6 









0 

C. 

1 

0 




1 













B 

1 

0» 

0« 

0-s 

Oft 

1 ft 

0* 

2ft 

1» 

Oft 

1ft 



C 

1 

0* 

04i . 

. Oft 

0* ■ 

1 ft 

Oft 

2 ft 

1ft 

Oft 

1 ft 



D 


0^ 

0^^ 

0« 

Oft 

1 ft 

Oft 

2ft 

Ift 

0« 

1 ft 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 


1 


6 

1 

0» 

0« 

0« 

0« 

1* 

C 

1 

0» 

0« 

0*. 

0* 


0 



o« 

0» 

Oa 

u 


B 

1 

0« 

Oa 

Qn 

0« 

1 « 

c 

1 

0«- 

0« 

Oa 

0* 

I * 

B 

1 

0* 

0» 

0» 

0« 

1 


1 


B 

1 

0* 

0« 

Oa 

o« 

1 «• 

c 


0* 

0* 

Oa 

0* 

la 


1 






c 


0* 

o». 

0-a 

Oa 

la 


1 


8 

1 

0* 

0« 

0« 

Oa 

1 a 

c 

1 

0* 

0« 

Oa 

Oa 

1 a 

0 

1 

Oa 

Oa 

0* 

Oa 

la 

0 


1* 

la 

"Oa 

l‘a’ 

■ 1 a 



LI 

L2 

L3 

LA 

L 5 

L6 


n 

3 

C 

0 



Oa 

3 a 

Oa 

Oa 

1 a 



Oa 

3 a 

Oa 

Oa 

1 a 



Oa 

3 a 

Oa 

Oa 

la 




LI 

L?. 

L3 

LA 

L ' 

L^- 


1 

0 

0 ■ 

2 



la 

Oa 

Oa 

2a 

la 



la 

Oa 

Oa 

2a 

1 a 



la 

Oa 

Oa 

2 a 

la 




LI . 

L2 

L3 

LA 

1 5 

L6 


1 

0 

1 

1 



la 

Oa 

la 

la 

la 



la 

0« 

la 

1 a 

1 a 



la 

Oa 

la 

la 

la 




LI 

L2 

L3 

LA 

L3 

L 6 


1 

1 

0 

1 



1 a 

1 a 

Oa 

la 

la 



la 

la 

Oa 

la 

la 



la 

la 

Oa 

la 

1 a 



Oa. 

Oa 

Oa 

Oa 

la 
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Table 

2-2. 

SERIES 

FOR RADIUS 

CONDITION 

(continued) 









ITI 

L2 ■ 

L3 

LA 

L5 

L6 








1 

1 

1 

C 




i 












8 

0* 

1' 

Oa 

0* 

Oa 

K 

la 

la 

1 a 

Oa 

la 



C 

1 

■0« 

Oa 

Oa 

Oa 

la 

la 

la 


Oa 

1 a 




1 












0 

0« 

Oa 

Oa 

0* 

la • 

la 

la 

la 

Oa 

la 




1 












0 


la 

la 

Oa 

la 

Oa 

Oa 

Oa 

Oa 

1 a 










LI 

L2 

L3 

LA 

L<3 

Lc 








2 

0 

0 

1 




1 












B 

Oa 

Oa 

Oa 

Oa 

la 

2a 

Oa 

Oa' 

la 

la 




1 












C 

Ok 

Oa 

o» 

Oa 

la 

2a 

Oa 

Oa 

la 

1« 




1 












D 

Oa 

Oa 

Oa 

Oa 

la 

2a 

0 a 

Oa 

la 

la 










Ui 

1-2 

L3 

LA 

L5 









2 _ 

0 

1 

0 




1 












U 

Z« 

Oa 

1 a 

Oa 

la 

Oa 

Oa 

Oa 

Oa 

la 










LI 

L2 

L3 

LA 

1 

L6 








■a 

0 

0 

0 




1 . 












B 

Oa 

6a 

Oa 

Oa 

J a 

3a 

Oa 

Oa 

Oa 

la 




1 












C 

Oa 

Oa 

0« 

Oa 

la . 

3a 

Oa 

Oa 

Oa 

1 a 




1 












0 

0« 

■ Oa 

Oa 

0* 

1 a 

3a 

Oa 

Oa 

Oa 

la 










U1 

L2 

L3 

L'‘ 

L " • 

L f- 








0 

1 

1 

2 




1 












c 

Oa 

Oa 

Oa 

La 

] a 

Oa 

1 a 

la 

la 

1» 










Ll 

L2 

L3 

LA 

15 

L6 








0 

1 

2 

1 




1 












0 

Oa 

Oa 

Oa 

la 

1 a 

Oa 

la 

2a 

Oa 

la 




1 












D 

Oa 

Oa 

la 

Oa 

la 

Oa 

1 a 

1» 

la 

la 
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Table 

2-2. 

SERIES FOR RADIUS 

CONDITION 

(continued) 










L 1 

L'2 

L3 

L^ 

Li; 

Li 









0 

X 

3 

0 , 




1 













0 


0« 

0» 

1«-- 

On 

1ft 

Oft 

Xft 

2ft 

Oft 

1ft 











LI 

L? 

L3 

xa 

L5 

. >-x- 









0 

2 

1 

I 




X 













c 

1 

0» 


Or 

Oft 

1 ft 

Oft 

Xft 

X» 

Xft 

1ft 



D 

1 

0* 

0« 

Oft 

In 

1 ft 

0« 

2 ft 


Oft 

Xft 



C 


0« 


0* 

Oft 

1 ft 

Oft 

X ft 

Xft 

Xft 

1ft 











LX 

L2 

L3 

L4 

IT- 

Lf 









0 

2 

2 

0 




1 













C 


0* 


0® 

Oft^ 

1 ft 

o« 

Xft 

2ft 

Oft 

Xft 




1 













c 

X 

0* 

. 0« 

I« 

On 

1 ft 

Oft 

2 ft 

Xft 

Oft 

Xft 



D 


o« 

X « 

0*. 

0 ft 

i ft 

Oft 

Xft 

2ft 

Oft 

In 











LX 

L2 

L3 

LA 

L-> 










0 

3 

0 ■ 

1 




X 













0 

X 

0« 

0« 

0-ft 

In 

1 ft 

Oft 

3 ft 

Oft 

0« 

1 ft 



0 


0* 

2 * 

0« 

Xft 

1 ft 

0« 

X ft 

0» 

Oft 

Xft 











L'l 

L2 

L3 

LA 

l5 

U< 









p 

3 

1 

0 




1 













c 

X 

0« 

x« 

On 

On- 

1ft 

Oft 

2ft 

1ft 

Oft 

Xft 



c 

X 

0< 


1* 

0ft . 

1 ft- 

Oft 

3 ft 

Oft 

Oft 

1« 



0 


0« 

x« 

0« 

Oft 

1ft 

Oft 

2 ft 

Xft 

Oft 

1 ft 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 










L’l 

L2 

L3 

L^ 

L5 ■ 

L6 









• 0 

4 

0 

0 




1 













c 

1 

0* 

U 

0* 

Ok- 

' IK 

0« 

3 k 

Ok 

Ok 

l.K 



D 


o« 


0 «• 

Ok 

1 K 

Ok 

5k 

Ok 

0« 

.IK 











LI 

L-1 

L-3 '■ 

L<1 

L 

L6 









1 

0 

0 

3 




1 









" 

' 



D 


o« 

0* 

0* 

IK 

1 K 

Ik 

._0« 

Ok 

■ 2k 

IK 






* 





LI . 

L2 

L3 

L4 

L‘='- 

L6 









1 

0 

1 

2 




1 











■ 


D 

.1 

0* 

0« 

Ok; 

■ iK 

1 K 

1« 

OK 

] K 

iK 

Lk 



D 


0* 

Oe 

1* 

OK 

i K 

Ik 

Ok 

Ok 

2 K 

Ik 











bj 

L2 

L? 

L4 

I'd 

L6 









1 

0 

2 

_i 




1 











- 


D 


Oit 

Ok 

Ik 

■ 0-K 

I K 

IK 

Ok 

. i 

Ik 

1 k 











L'l 

L? ■ 

L3 

L4 

15 - 

L* 









1 

1 

0 ' 


‘ 



1 













C 

1 

b* 

Ik 

Ok 

Ok 

1 * 


Ok 

■ "b K 

“^K 

IK 



D 

•1 


Ok 

Ok 

Ik 

1* 

Ik- 

Ik 

"Ok 

'• Ik 

iK 



D 


b* 

1 « 

0 n 

- 

1 K 

Ik 

Ok 

Ok 

2k 

l« 











U1 

l2 

L3 

L4 

Lf 

it' 






, 



1 . 

1 

1 

1 ' 




1 













c 


0* 

1« 

Oc 

Ok 

1 K 

‘1» , 

■ Ok 

Ik 

IK 

IK 




1 













C 



0» 

0« 

0« 

iK. 

Ok 

Ik 

Ik 

IK 

IK 




1 













D 


o« 

0« 

Ok 

1 » 

1 * 

Ik 

1 

1 K 

Ok 

i« 




i 













D 

i 

0^ 

0* ' 

1 K 

0» 

1 K 

1« 

IK 

i 

Ok 

Ik 

IK 



D 


0* 

IK 

. 0* 

Ok 

1 K 

iK 

Ok 

1 K 

Ik 

Ik 




1 



. 

■ ■ ■ 









D 


i» . 

Ok- 

Ok 

Ok- 

1 K 

• Ok 

. IK 

iK 

1« 

IK 
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Table 

2-2, 

SERIES 

FOR RADIUS 

CONDITION (continued) 









L-l 

I2 

L3 

Lo 









I 

1 

? 

0 


1 











c 

1 

i« 

0« 

0» 

Os 

Is 

Os 

Is 

2s 

Os 

Is 

D 

1 

0* 

0« 

I< 

0* 

1« 

Is 

1 « 

1 s 

Os 

1 s 

D 


1* 

0< 

Os 

■ Os 

1« 

Os 

Is 

2« 

Os 

1 s 









El 

L2 

L3 

LA' l5 









1 

2 

0 

1 


1 











C 

1 

0» 


Ok 

ps 

1 s 

Is 

Is 

Os 

Is 

Is 

0 

1 



OS 

Os 

is 

Is 

Is 

Os 

Is 

1 s 

D 

1 

0* 

c « 

Os 

1* 

1 s 

u 

Os 

Os 

Os 

1 « 

D 



1 ■* 

Os 

L ^ 

1 s 

Os 

1 s 

Os 

Os 

Is 









LI 

L? 

L3 

La L3 










2 

1 

0 


1 











C 

1 

0« 


Os 

Os 

i * 

1 S 

Is 

is 

Os 

IS 

C 

1 

1# 

0* 

Os 

Os 

1 s 

()S 

2s 

is 

Os 

is 

c 

1 

0* 

1* 

Os 

Os 

1 s 

Is 

1 s 

is 

Os 

Is 

c 

1 


Q^^ 

Os 

0« 

i s 

Ok 

2 « 

1 « 

Os 

Is 

0 


1* 

1« 

Is 

Os 

Is 

Os 

1 s 

0* 

Os 

is 









LI 

L2 

L3 

LA L^ 









1 

3 

0 

0 


1 











c 

1 


0* 

Os 

Os 

Is- 

Os 

1 * 

Os 

Os 

is 

c 


1^. 

0» 

Os 

Os 

is 

0* 

^ s 

Os 

Os 

• 1 s 
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Table 2-2. . SERIES FOR RADIUS CONDITION (continued) 










LI 

L2 

L3 

lA 









2 

0 

r 

2 


1 











c 

1 

1* 

0« ' 

o» 

0« 

la 

la 

Oa 

Oa 

?a 

la 

c 

1 

0* 

o« . 

0« 

■ 

i « 

2 a 

Oa 

Oa 

la 

1 a 

c 


1* 

0* 

0« 

0< 

1 » 

1 a 

Oa 

Oa 

.2* 

la 









Ll 

L2 

L? 

LA 









2 

0 

1 

1 


1 











c 

1 

1* 

0«‘ 

0^^ 

0* 

la 

la 

Oa 

la 

la 

la 

D 

1 

0^ 

0« 

U ■ 

0« 

1 a 

?* 

6a 

Oa 

la 

1 a 

C 


1* 

0« 

. Ofi 

0» 

I a 

1 a 

Oa 

la 

la 

1 « 









U1 

L2 

L3 

LA 









2 

1 

0 

1 


1 











C 

1 


I « 

0* 

0* 

!•» 

2a 

oa 

Oa 

1 a 

1 a 

C 

1 

1* 

■ 0* 

0-H; 

Oa 

1 * 

1 a 

1 a 

Oa 

la 

la 

D 

1 

0* 

Hf 

Q< 

Oa 

I a 

2a 

Oa 

Oa 

la 

1 a 

0 

1 

1» 

c« 

Ot 

0^^ 

i a- 

la 

1 a 

Oa 

la 

la 

D 


1 « 


0^ 

1-S 

I* 

La 

Oa 

Oa 

Oa 

1 a 

• 








LI 

L2 

L3 

*-A 









2 

1 

1 

0 


1 











c 



0* 

0» 

0-2 

i * 

la 

1 a 

la 

Oa 

1 a 


1 











rs 

ly 

L 

l-« 



0^ 

1 * 

la 

1 a 

la 

Oa 

1 a 

D 


♦ 

1* 

la 

0« 


la 

Oa 

Oa 

Oa 

1 a 

D 


z* 

0« 

. 


la 

Oa 

1 a 

Oa 

Oa 

la 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 









LI 

L2 

L3 

- L4 L5 


i 







0 

0 

1 

c 

i 

0* 

Oft 

Oft 

1ft 

C. ft 

Oft 

Oft 

ift 

1 ft 

c 

Oft 

1 

Oft 

Oft 

1 ft 

Ift 

?« 

Oft 

Oft 

Oi 

1 ft 

D 

1ft 

' Oft 

Oft 

Oft 

1ft 

2 ft 

0 ft 

0 ft 

1 ft 

1 ft 








LI 

L2 

L3 

LA ir 


1 







0 

1 

0 

D 

0 ft 

l 

0 ft 

la 

Oft 

ift 

1ft 

Oft 

Oft 

Oft 

Ift 

D 

2* 

Cft 

1* 

Oft 

Ift 

Ift 

Oft 

Oft 

o» 

1ft 








Ll 

L? 

L3 

L4 L‘> 


1 






■j 

1 

0 

0 

C 

Oft 

1 

Ift 

Oft 

Oft 

1 ft 

1ft 

Oft 

Oft 

Oft 

1 ft 

D 

Oft 

.1 ft 

Oft 

0 ft 

1ft 

3ft 

Oft 

Oft 

Oft 

1 ft 


, 






Li 

L2' 

L3 

L4 


1 




- 



0 

0 

0 

C 

1» 

1 

Oft 

0* 

Oft 

1 ft 

3ft 

Oft 

Oft 

Oft 

1 ft 

D 

1ft 

Oft 

0ft 

Oft 

•i * . 

3ft 

Oft 

0* 

Oft 

1ft 





... 


11 L2- I 

-3 L4 l5 


1 






i 

2 

0 

2 

B 

Oft 

1 

Oft 

0-ft 

Oft' 

I ft 

1ft 

2 ft 

0* 

2ft 

Ift 

C 

Oft 

1 

Oft 

Oft ' 

Oft 

1 ft' 

Ift' 

2ft 

O'ft' 

2 ft' 

1«‘ 

D 

Oft 

0 « 

Oft" 

Oft 

1» 

1ft 

2 ft 

Oft 

2ft 

1ft 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 

Ei L2 L3 L'' L6 

1 ■ 1 1 


B 

0* 

1 

0* 

Os- 

Os 

1 S 

Is 

2 s 

Is 

Is 

Is 

C 

0* 

1 

o« 

Os 

Os 

i * 

1« 

2 s 

1 S 

Is 

1 S 

D 

0^ 

0^ 

Os 

Os 

I * 

Is 

2 s 

is 

.Is 

is 








LL 

L? 

13 

LA 


1 






1 

2 

2 

0 

8 

0» 

1 

0« 

Os 

Ok 

i » 

is 

2s 

2s 

OS 

Is 

C 

0* 

1 


o« 

Os 

i * 


2s 

2s 

Os 

Is 

D 

0« 

0« 

0* 

Os 

Is 

Is 

2s 

2s 

Os 

is 








U1 

L2 

L3 

L4 L = 


1 






c 

1 

0 

2 

B 

0« 

1 

0^ 

Ds 

Os 

1* 

2s 

Is 

Os 

2» 

Is 

C 

0* 

1 

0« 

Os 

0« 

is 

2« 

is 

Os 

2« 

1 s 

D 

0» 

o.» 

Os 

0* 

. 1« 

2 s- 

Is 

0 s 

2« 

Is 








L'l 

L2 

L3 

LA L5 


1 






2 

1 

1 

I 

B 

0«- ■ 

1 

0* 

0« 

Ok 

i » 

2S 

Is 

1 s 

1 s 

Is 

C 

0« 

1 

0« 

0« 

Os 

]S 

2s 

1 s 

Is 

Is 

Is 

D 

0«- 

0*. 

0* 

o» 

Is 

2s 

Is 

Is 

Is 

Is 







Ll Lz Ls la L5 


1 


- 




2 

1 

2 

0 

B 

0* 

I 

0* 

Os 

Os 


2s 

1 s 

2s 

Os 

1 s 

C 

o» 

1 

0 ^ 

Os 

o» 

is 

?s 

is 

2s 

Os 

Is 

D 

0* 

0* 

Ok 

0* 

Is 

2 s 

1 S 

2« 

Os 

1 S 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 










L'l 

L2 

L3 

L4 

1. 3 

Lf 









C 

3 

1 

2 




1 













c 


0« 

2« 

0« 


1 * 

0« 

1 « 

1* 

1« 

1 « 











LI 

L2 

L3 

L4 

L3 

K 









0 

3 

2 

1 




1 













D 


0» 

2< 

o« 


1 * 

0« 


Z* 

0« 

1« 











LI 

L2 

L3 

L4 

L5 

L6 









. 0 

4 

1 

1 




1 













D 


0* 

1 « 



1* 


2 * 


0» 

1« 











Ll 

L2 

L3 

LA 

Lf 

L6 









0 

5 

0 

1 




1 













D 



2« 

0* 

L * . 

1 


3 « 

0<t 

0« 

1 -K 











U1 

L2 

L3 

LA 

Lf 

L6 









1 

2 

0 

3 




1 













D 

1 


d« 

0» 


] <- 

1« 

2« 

0« 

2« 

1^* 



D 


0*>- 

2 « 

0* 

1* 

1 «■ 


0* 

0* 

2« 












Ll 

L2“ 

L2 

LA 

L5 

L6 









1 

2 

1 

2 




1 













D 

1 

0* 

0» 

0« 


l^i 

1* 

2 « 


m 

1 « 



D 

1 

0« 

0« 

i« 

Oe 

1* 

l«‘ 

2 e 

0* 

2« 




0 ■ 



2« 

0« 

1 

1* 

1* 

0« 

1« 






1 

' 












D 



1« 

0» 


1^* 

0« 


1* 

!«■ 

1 « 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 










Li 

L2 

L3 

L^ 


L6 









1 

2 

2 

1 




1 













c 

1 

Ofr 

’ 0« 

Oa-' 

I'a 

i'a 

la 

?. a 

2 a 

Oa 

la 



D 

1 

0« 

G « 

u 

Q-a 

la 

la 

2 a 

1« 

la 

la 



C 

1 

1« 

■ 1«. 

Oa' 

1* 

1 a 

Oa 

1 a 

2« 

Oa 

la 



c 


1 » 

1 6 

1 » 

Qa 

1 a 

Oa 

1 a 

la 

. 

1 a 











LI 

L2 

L'3 

1.2 

1 =: 

L^- 









1 

2 

3 - 

0 




1 













D 

1 

0^^ 

0« 

la 

Oa 

1 a 

1 a 

2 a 

Z a 

Oa 

1 a 



C 


1* 

Is 

1 * 

pa 

la 

Oa 

la 

2a 

Oa 

1 a 











LI 

L2 

L3 

L2 

1“" 

Lfc. 









1 

3 

b 

2 




1 













c 

1 

0*^ 


Oa 

0* 

1 a 

la 

2 a 

0* 

2* 

la 



c 

1 

0« 

T« 

0« 

Oa 

la 

la 

2 a 

Oa 

2a 

la 



c 


0^ 

2 a 

Oa 

La 

1 a 

1 a 

la 

Oa 

la 

la 











LI 

L? 

L3 

LA 

L*^ 

Lf 








. 

1 

3 

1 

1 




1 













c 

1 

0* 

la 

0« 

0* 

1 a 

la 

2 a 

1 a 

la 

la 



p 


0* 

la 

0* 

Ot 

1 a. 

1 a 

?’a 

1 a 

la 

la 




1 













c 

1 

0« 


0« 

la 

I a 

la 

1 a 

1 a 

0» 

1 a 



0 



la 

o« 

la 

1 a 

Oa 

2 a 

1* 

Oa 

1 a 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 










Ll 

. 

L3 

L4 

L3 

L6 









1 

' '3' '■ 

■”2 

0 




I 













c 

1 

0* 

1« 

Ok 

o» 

1« 

IK 

2k 

2k 

0« 

1 « 



D 

I 


1« 

0* 

Ok 

1 K 

iK 

2 K 

2 K 

0« 

1« 



D 




■ 

o« 

•IK 

0« 

2# 

!«■ 

Ok 

IK 











Ll 

L2 

L3 

L 

L£ 

L6 









1 

A 

0 

1 




1 













D 


1* 


Ok 

1* 

3 * 

o« 

3 K 

o« 

Ok 

iK 











Li 

L2 

L3 

L4 

L5 

L6 









1 

4 

1 

0 




1 













D 



u 

1 « 

&K 

] » 

Ok 

-3 « 

0* 

C« 

iK 











Ll 

L2 

L3 

L4 

L'j 

L^’ 









2 

1 

G 

3 




1 













C 

1 



Ok 

I Ji- 

1* 

2k 

IK 

Ok 

2k 

1» 



D 


1* 


Ok 


i » 

1« 

Ok 

Ok 

2k 

i K 











Ll 

L2 

L3 

L4 

13 

L6 









■V 

c 

1 

1 . 

2 




1 













0 

1 


0* 

Ok 

IK 

IK 

2« 

iK 

1 » 

Ik 

l« 



D 

1 

o« 

Ok 

l< 

Ok 

Ik 

2« 

1 « 

0« 

2^. 




C 

1 

1* 

i« 

Ok 

i« 

I K 

Ik 

o« 

l« 

l« 

IK 



D 


1« 

1« 

1-K 

Ok 

1 * 

Ik 

0* 

0« 

2« 

1* 




'2-28 


TR-792-8-283 


NORTRONICS - HUNTSVILLE 




Table 

2-2. 

. SERIES 

FOR RADIUS CONDITION (continued) 









LI 

L2 

L3 

LA 15 









2 

1 

2 

1 


1 

- 










D 

i 

Q* 


0« 

Is 

Is 

2S 

1 S 

'2s 

Os 

is 

D 

1 

o« 

0« 

1* 

. Os 

is 

2 * 

is 

1 s 

Is 

1* 

C 


1» 


u 

Os 

1 * 

is 

Cs 

is 

is 

is 


1 




- 







D 



0« 

1 « 


Is 

Os 

is 

1 s 

Is 

Is 






■ 



LI 

L2 

L3 

LA 15 









2 

1 

■3 

0 


1 

i 










0 

1 

0» 

0« 

1« 

0« 

Is 

2s 

Is 

2s 

Os 

1 s 

0 


2 


1« 

Os 

! s 

Os 

is 

2s 

Os 

1« 




- 


- 



LI 

L2 

L3 

LA 15 









2 

2 

0 

2 


1 











c 

1 

0« 


Os 

Os 

1 S 

2s 

1 » 

Os 

?« 

Is 

c 

1 


0» 

Os 

o» 

is 

Is 

2 s 

Os 

2s 

is 

D 

1 

0* 

is 

Os- 


1 S 

2s 

is 

Os 

2s 

is 

0 

1 

0* 

2 * 

Os 


is 

2s 

Os 

Os 

Is 

1 s 

c 

1 

1« 

C« 

0 * 

Os 

is 

Is 

2« 

OS 

2s 

is 

0 


1* 

1 « 

Os 

is 

!s 

Is 

1 s 

Os 

is 

is 









LI 

L2 

L3 

LA L5 









2 

2 

1 

3 


1 











c 

1 

0* 


0« 

0* 

l‘s 

2s 

is 

Is 

1« 

1 K 

c 

1 

l» 

0« 

Os 

Os 

is 

Is 

2« 

Is 

Is 

1 « 

8 - 

I 

0» 

i« 

Os 

0's‘ 

1 S 

2« 

is 

is 

is 

is 

C 

1 


0-k" 

Os 

Os' 

I's 

is 

2s 

is 

is 

is 

0 

1 


1« 

Os 

1 S 

Is 

Is 

is 

■ is 

Cs 

1» 

D 



1 » 

Is 

Os 

1 s 

Is 

is 

Os 

1 « 

1 s 
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Table 

2-2. 

SERIES FOR RADIUS 

CONDITION 

(continued) 










LI 

L2- 

1.3 

■ 

L? 

L6 









2 

2 

2 

0 




1 













c 

1 

0* 

1 « 

0« 

Oft 

1 ft 

2ft 

1 ft 

2ft 

Oft 

1 ft 



c 

1 

1» 

0« 

o» 

Oft 

1 ft 

1ft 

2 ft. 

2ft 

Oft 

1ft 



0 

1 

0* 


0* 

Oft 

1ft 


1ft 

2ft 

Oft 

1 ft 



D 

1 

1 * 

0« 

Oft 

Oft 

1ft 

1ft. 

2 ft 

2ft 

Oft 

Ift 



0 

1 


n* 

L* 

Oft 

1ft 

1ft 

1ft 

1ft 

Oft 

1ft 



D 


3* 

0« 

Lft 

0< 

1ft 

Oft 

2« 

1 ft 

Oft 

1ft 











LI 

L2 

L3 

LA 

L5 

L& 


T 







2 

3 

1 

0 



!3 

1. 

2 « 

0« 

1ft 

Oft 

Ift 

Oft 

3ft 

Oft 

0« 

1 ft 











LI 

L2 

L3 

LA 

L5 

Lf- 









"! 

0 

1 

? 




1 













D 


3# 

0« 

i ft 

0 ft 

1ft 

1ft 

Oft 

Oft 

2ft 

1 ft 











Li 

L2 . 

L3 

LA 

L5 

L^, 









•5 

0 

2 

1 




1 













D 


2* 

0« 

U 

0ft 

1 ft 

1ft 

Oft 

1 ft 

1ft 

1 ft 











LI 

L2 

L3 

LA 

L 5 

.L^ 










1 

0 

2 




1 













C 

1 

1* 


0^ 

Oft 

1 ft 

2ft 

1 ft 

Oft 

2ft 

1 ft 



D 

1 

1* 

0« 

Oft 

Oft 

1 ft 

?ft 

1-ft 

Oft 

2ft 

1 ft 



D 


1» 


Oft 

1ft 

Ift 

z« 

O-ft 

Oft 

1« 

1ft 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 


LI L2 L3 LA L6 

• 3 1 I 1 

1 


c 

1 

1« 



Ofi 


D 

1 


0“ 


0« 


D 

1 

l» 

I 

u 

0* 

1 * 

D 


Z* 

0^^ 

I* 

0» 

J* 


c 

1 

1* 

0« 

0» 

0® 

i a 

D 

1 


qa 

oa 

Oa 

la 

0 

1 

2a 

0* 

la 

0» 

la 



1 






0 

■ 

0* 

2 * 

Oa 

la 

1 a 


1 






c 


2a 

0* 

1 B 

Oa 

la 


L 






D 



la 

Oa 

la 

la 


1 






0 


la 

L'a 

la 

0« 

la 


?a 

1 a 

I a 

1 a 

la 



2» 

la 

la 

la 

la 



2 a 

Oa 

0 a 

la 

l« 



la 

1 a 

Oa 

la 

la 




LI 

L? 

Ll 

LA 

L5 

L{- 


3 

1 

2 

0 



2a 

la 

2a 

Oa 

la 



2 a 


2a 

Oa 

la 



la 

la 

la 

Oa 

la 




111 

L2 

L3 

L^-* 

L5 

L6 


3 

2 

0 

■1 



3 a 

Oa 

Oa 

Oa 

l.« 




11 

L2 

L3 ■ 

•LA 

LF' 



A 

0 

1 

1 



2 a 

Oa 

Oa 

la 

la 




U1 

L2 

L3 

LA 

Lp 

L6 


A 

1 

0 

1 



3 a 

Oa 

Oa 

Oa 

1 a 




U1 

L2 

L2 

LA 


Ld- 


A 

1 

1 

0 



3a 

Oa 

Oa 

Oa 

la 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 









Ul 

L2 

L 3 

14 

15 

L6 


l' 






- 

0 

1 

0 



c 


0« 

I'ft 


i# 

3» 

0-« 

O'#- 

0« 

1# 










Ul 

L2 

■ L3 

L4 

LT 

L6 


1 






1 

4 

0 

3 



c 

0* 


0* 

1 «■ 

u 

1* 

2« 

0» 

2» 

1 










Ul 

L2 

13 


L- 

L4 


1 






1 

4 

1 

n 



D 

o# 


0* 


1 » 

u 

2# 

1* 

1« 

1 « 










Ul 

L2 

L3 

L4 

L5 

ir. 


1 






1 


2 

1 



0 

0* 

2* 

o« 

1 « 

i» 


2 » 

2« 

0« 

1 ■« 










Ul 

L2 

L3 

14 

f_5 

It 


1 






2 

3 

0 

5 



c 

1 

2-m 



I 

2* 

1 « 

0« 

2« 

1^* 



D 


i-« 

0» 


1* 

1« 

2# 

0# 

2« 

1« 










Ul 

L2 

E'i 

14 

14 

Lt 


1 






2 

3 

1 

2 



0 

. -O'* 

1 

2 « 

o« 


1 * 

?» 

1 « 

1* 

1« 

1« 



D 

1« 

I 

1« 

Ofc 

1-if. 


1# 

2# 

1# 

1« 

l^^ 



D 

1» 

If 

1« 

O'# 

1* 

1« 

2# 

0« 

2# 

L« 

• 
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Table 2-2. SERIES FOR RADIUS CONDITION (continued) 










L'l 

L2 

L3 

LA 

L 5 

L6 









tl 

3 

? 

1 




1 













D 


0* 

2*- 

0< 


1 

2« 

1 « 

2* 

04* 

144 




1 





' 








C 

1 

1» 


0* 

1« 

1* 

1« ' 

2^t 

2» 

044 

144 



D 


1* 

1« 

■ 1* 

0* 

1 * 

u 

2^< 

i* 

14! 

144 











LI 

L2 

L3 

LA 

L3 

L6 









2 

3 

3 

0 




1 













D 



1« 


0» 

1 « 


2« 

2» 

044 












LI 

L2 

L3 

LA 

L5 

L( 









a 

2 

0 

3 




1 













D 


1* 

1 « 

0*- 

i« 


2 « 

1* 

0» 

2 44- 

144 











LI 

L2 

L3 

LA 

L3 

Lf 









3 

2 

I 

2 




1 













0 




Os 

1* 


2^> 

i * 

1* 

l4l 

l4i 




1. 













C 

1 

' 1* 


1« 

0-^ ’ 

1 ^ 

2« 

l« 

04f 


l44 



D 


2* 


_ 

a* 

1 


2 « 

0* 

244 

l4i 











LI 

L2 

L3 

LA 

1 5 

L<^ 









3 

- 

2 

1 




1 













R 

1 

1 » 

1* 

0« ' 

i« 

1 

2» 


2* 

044 

144 



D 

1 

1 * 

1« 


o« 

1 « 

?. « 

1 « 

!«■ 

■ 1 4! 

I 41 



0 


2« 

0« 

’ .1 « 

0* 

1 «• 

1» 

2» 

1 * 

144 

l44 











Ll 

^2 

L3 

LA 

L? 

L6 









3 

2 

3 

0 




1 













D 

1 



1* 

0* 


2* 

i« 

2 4^ 

044 

1 44 



C 



0« 

1* 

0» 


1« 

241 

2» 

044 

l4t 
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Table 2-2. SERIES FOR RADIUS CONDITION (concluded) 


LI L2 L3 L4 L5 L6 
^ 11 2 


D 

2» 

Q* 


0« 

I » 

2^ 

i If 


2 « 

l^i 





- 




LI 

L2 

L3 

L4 L = 

L6 








c. 

1 

2 

1 



I 











C 

■ 2« 

0« 

1.*-. 

OV 

I « 

2* 

1« 

1* 

1» 

1» 









LI 

L? 

L3 

L4 L r 

L6 








4 

1 

2 

0 



1 











D 


0« 

1* 

0<- 

1 R 

2 * 

1 « 

2 » 

0« 

1« 



SAVE TAI-'c 

CASE lu= f %ir-2-RC 0*«2 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION 






— 

— 

— 

-Irl-— ■ 

Trz — 

~tri — 

'•t'M'" ' 





. 



1 

C 

C 

c 

A 

X ■ ■■ ■ 

C» 

0» 

0« 


L«. 






B 

1 

c ■** 

c« 

c « 

c* 

1* 

o« 


c* 

c « 


C 

T 

C* 

c« 



1* 

o» 

.]« 

c « 

0 « 

1« 

c 

1' 

c« 

c« 

r * 

C« 

1« 

0» 


c « 

c « 

1« 








L 1 

l2 

L3 

L4 









. T. . 

-- 

A 













1 










H 



r-* 


b* 

0* 

!-«- 


G« 



1 










B 

€ * 


-j* - 



0« 

-3« 

C^> 

C- « 



1 


















r* If 


— t-H- 

6“ 

— 



' L 

1. « 

L' * 





1 










- --- c 


— . 

0-^ - 

-G-*- 

1*- 

0« 

• C-« 

• c 

• 



1 










D- 

■ fy-Vr 

- 

'0^'* 


- !•«• 




G^» 







1 K 

p tf 



C » 


B 


i-R — 



i *♦ 













L 1 

L2 

L3 

LA 








1 

1 

r 

C 

C 


1 










B 

c « 

c« 


C* 

1« 

. 

c ^ 

L « 

C « 

1« 

B 

1 » 

c « 

c» 


m 

0 

C* 

C « 

C« 

1« 


L 










C 


0^« 

c« 

c« 

i« 

1 « 

0 * 

C^> 

c « 


C 

• 1- - ■ 
i« 

c •» 

0» 

0# 

i« 

0» 

0# 

c « 

0 « 

i« 


1 










c 


c « 

0 ^ 

c« 

i« 

1« 


c« 

c « 

i« 


1 -■ 

- ■ 

- 








c 

1 « 

Q « 

c« 

c« 


0^f 

•J« 

c« 

c 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 









Ll 

L2 

L3 

LA 


~ - — 


— ^ 

— 

.... . . 

, 

^ . 


... . 

n 


1 










e 

r. » 


0« 

C« 


0 » 

1« 

c « 

C R 

IR 


1 










• c ■ 



On 


i« 

0« 

i« 

c * 

■ c « 

IR 


1 










c 


— 

C 



-Oir- 

— 1-1! 

” C-»- 

-■e-R 

-IR- 


1 










c 

C » 

2 


c« 


C« 


c « 

o« 

i« 




- 




L r 

L2 

L3 

LA 








1 

1 

C 

r 

B 

c * 



c« 

1* 

l< 

3-» 

c » 

c « 

i« 


1- 










B 

1 » 




1^1 

3« 

1« 

c * 

c« 

iR 


1 


- 





- 



C 

'■ a 

1 * 

c 

:•» 

i« 

l« 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 


2 — ■ 't L'6 
1 1 1 


• 1 


B 

i 

C » 

c « 

C« 

C* 


2 * 

1* 

1« 

1 « 

1« 

C 


C« 


0* 

0 < 

1« 


1* 
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Table 2-3. SERIES FOR ORTHOGONALltV CONDITION (continued) 


Li L2 L3 LA Lf L6 
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B C« C« L» 0« 3* C« C« I-# 
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C C« 0« 0* 1^* 

I • 
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■ 1 

D C« 2« C* 1« 0« l-» C« Q« 1« 
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c- C* - 'i« . 0« l^s 0* O* C » C« 1« 

Li- -Lx- ■ t-3 -.LA • -ti L6 

1 C C 2 

i .. .. 

8 C* C* C« C* 1* 1« C« 2« l« 

C C» C« C-» C« 1« l^* L« C« 2 * 1« 

C C« C« C-» C« 1-R 1« C« C« 2» 1« 

C 1« f* C« 2< l-« S« C«‘C« 1» 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continiied) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 


LI L2 L2 LA L£ 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES. FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITIOH (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES F.OR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 










LI 

L2 

L3 

iH 


1 







■r— 

■ - 

u 

Z 

B 

1 

« 

L« 

Oa 

0« 

la 

2a 

ia 

C a 

<1 « 

la 

• B 

1 


-c* 

Ca 


i a 

la 

2 a- 

C a 

2 a 

la 

C 

1 

ClF ' 

•■■1 It - 

‘ G* 

C'a 

la 

2« 

la 

C a 

2 a 

ia 


1 

1 •** 

C a 

0 ^ 

■“Ga — 

—ta- 

ira— 

— 2a 

— G a— 

-- 2-a- 

— 1«- 

C 

1 

1» 

I« 

n« 


• la 

la 

la 

Ca 

1 a 

la 

D 

1 

t: a 

1« 

0* 

•Oa 


2a 

1* 

C a 

2a 

la 


1 

IL » 

L * 

~C a' 

-G* 

1T( 

ta— 

-2a-- 

— G"a — 

■■■2a - 

-l-« • 

D 

1 

1«- 


o«- 

la 

la 

la 

■la 

C a 

ia 

la 

D 

p — 

1 

1 

-1-5 

2a- 

-0 -B. 

Ga - 

-ea- - 


la 

la 

■T-« 

• G« 

2 a 

la 


1 

-C - - i-sj .Q if 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 



’ " — 

— - — 

• 

- - . 

— - 


tl 

L2 

L3- 

LA 


. — 






2 

2 

1 

1 

B 

1 


C-s 

C« 

l« 

2« 

1« 

i« 

la 

la 


1' 

- ' - 

- 

■ 







B 


c« 

C« 

c* 



2* 


la 

- la 


1 ^ 



• '■* 




... 



C 

c« 

1* 

0« 

c* 

U( 

2« 

1-*^ 


la 


C 


2 « 

1* 

c« 


2« 

C« 

c« 

1 a 

la 


1 

• - 


- 







C 

1 

C« 

c« 

c« 


1 * 

2« 

1 « 

1 * 

la 


X 










C 

1 » 

— — — 


o« 


m 

1* 

1* 

1 « 

C« 

la 

C 


1* 

L« 

r« 

i« 

1« 

1» 

c « 

la 

1« 


i ■ 



‘ 




- 

. 


c 

2* 

c-» 

c» 


ifl 

0« 

2» 

1 « 

Ca 

la 

- . . 

1-" 


•“ - 

— - 

- 






c 

c-» 

— 

l« 

c« 


1* 

2« 

1« 

. u* 

la 

la 

D 

jCL 

c * 

2« 

1« 

c« 

i* 

2« 

o« 

c « 

la 

la 


. - 1 — 

■ 

— 







- 



c 

1 « 

c * 

0« 

c» 




i « 

la 

la 


~ - 



- 


. ^ - 





D 

‘ 1« 






i» 


Ca 

la 

C 

1« 

1 » 

1« 

c< 


1« 


c« 

la 

la 


1 

• - 




' 





c 

1« 

2« 

0« 

0« 

1* 

1* 

0* 

1« 

i a 

la 


1 


- 








D 

—5 

2« 

l^» 

1« 

1* 

1« 


c • 

Ca 

la 

0 


C « 




0« 

2» 

1 a 

Oa 

1« 

- — 

2- 



. 

. 






c 



G« 


1* 

0« 

1» 

i « 

la 

la 


-i 

. 

— 

~ * 







c 

2* 

!•» 


1* 


0* 

1« 

c « 

Oa 

la 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 










LI 

L2 

L3 

L^ 




~ 

• 

- 




i. 

2 

C 


1 











B 

1 

c-« 

i« 

■ 0* 

■G« 

1* ■ 


1* 

2 * 

. c« 


B 

1 

1 

<z « 

L « 

C * 

1« 


i * 

2 ^ 

0« 

1« 

C 

1 

■i: » 

!•* • 

'^J » 

0* • 

1* 

?.« 

• u 

' z-« 

0 * 

1* 

C 



- ' 



1« 


2^^ 

2 * 

c « 

1« 


-lr« — — 


-t « I"* — &7i» — i. 


C’ « 1 «■ 

0« 

0^ 

1« 

2« 


2 * 

0-» 


1 « c~» 

C-« ■- 

c-« 

i« 

I'S 

2» 

2 « 

c« 

< --M 

TTf TTf” 


-tjTF- 

1-5- - 

-1-Tf- 

^1-5“ 

1 «- 

C^5 - 

1 « 

1-^t — 

-2^— 


- 1«- 

- l^r- 


-C-5 

'C^• 

1« 



- -^TT - 


— i-5- 

0 5' 

1* 

2 « 


i « 


)T 


LI L2 L3 L4 L5 L6 
2 3 C 1 


c 

2» 

0 « On 

In 

i« 

0 ^ 

3 * 

C « 

Cn 

in 

D 

c ■» 

35 c« 

0 ^ 

m 

2n 

On 

Cn 

In 

in 

C 

1 

In 

2n Cn 

0« 

i« 

m 

in 

0 n 

1 n 

in 

C 

1 

2* 

C n 0 « 

In 

m 

On 

jn 

Cn 

On 

i « 

ir- 

L 

1-^— 

-2Tr-*“ay ■ 

-Ctt- 

-Ln 

1-5 

LI 

2 

1« 

L2 

In 

L3 
- 1 

Cn 

L^ 

C 

in 


I- - l-v ~ O--# 1 « '■ - • 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 





... . — 

— 

. 

.. - 

- 

LT 

L 2 -L3 









2 

A C 

t-' 

•v 

« ~ 

2~ 

.. 

— 

— 



... 



- 


D 


1* 

C« 

Oa 

1 a 

Da 

3* 

C a C a 

1« 








Ll 

L2 L3 

LA 


- - 

. - 

— 


- 


3 

c e 

■3 


1 









- • c - 

- 2^- 

tj< ■- 

tjTf 

■ I* 

la 

la 

t‘a 

• C a 2 a 

ia 


1 







— 

— 

CT' 

— rr— 

— 

L-ft 

2« 

—fa 

■'2 a 


^ L ^ 

4 - ^ 

C" - 

1 

2’T- 

t a- 

Ga 

-ta ■ 

la 

la 

■Ca 

e a 2 a 

la 



.. 





Ll 

•L2 L3 









3 

C 1 

2 

c 

"t 


Ga 

La 

la 

la 

C a 

la la 

la 


1 









c 

1« 

c« 

la 

la 

la 

2a 

■ V a 

C a la 

ia 

- 

1 - - 









c 

2 » 

Ca 

Qa 

la 

la 

la 

ca 

la la 

l« 


Ll L2 L3 LA II 
3 1 C 2 


• ” B ■ 

L 

1-a 

C a-- 

Ca- 

■G-a 

la 

2* 

•la 

C-a 

2a 

la 

C 

1 

1 a 

Ca 

Ca 

G» 

i» 

2 a 

1* 

ca 

2a 

la 

c 

1 

1 a 

1 a 

Oa 

,la 

la 

2a 

C« 

Ca 

la 

la 


1 











- c 

1 

2 a 

C a 

Ca 

la 

la 

-■T 

1 a 

la 

Ca 

-{i ,f--- 

1 a 

Tjf- 

la 

-la 

Cf 


-r«~ 

C * 

— GT5~ 

C a 



l * 


c. 

0 

1 

1 a 

1 a 

Ca 

1-a 

i « 

2« 

Qa 

C a 

la 

la 


I 











-- D- 


2" a- 

-C-a 

Ca 

-ta 

la 

la 

l« 

C a 

la 

l« 


2 




-*■ 


— T— * 

— — 

-C"«— 

-2-a 

-I« 

tT" 



A "a 

C-a 

C a i •“ 



C. " 


1 











C 


?a 

la- 

Oa 

2» 

la 

la 

Ca 

C a 

C a 

■ l« 


L6 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 



_ „ 

. 

. — 

— 



LT 

Lz 

Lt3 

LA L 








- 

1 

1 

1 

B 

■ 1« 

•j * 

G« 

C« 

u 

2* 


1 « 

- 1« 

IM 


1 • • 

-- 

" " 








C 

1* 

c« 

c « 

G« 

l« 

2* 

1* 


1« 

IM 

C 

1 

1« 

L« 


C« 

u 

2« 

C« 

C« 

1« 

IM 

C 

— t 

0» 

G* 


1* 

1* 

1* 


C« 

IM 


1 










C 

C* 

1 « 


1« 


3« 

U * 

G« 

c « 

iM 

c 

1 ■■ " 

0« 

0» 

c« 

l« 

2« 



1 «■ 

1« 

c 


1* 

1* 

0« 

1^* 

2« 

G* 

C « 


IM 



■■ 1- 

— 



- — ~ 


- 

- “ ■ 

* 



D 

2« 

c * 

ri-ti 


1« 

1« 


1* 

c « 

1« 



2- 




— 

" - 


- 




D 

2« 

1^» 



1« 

1« 

0« 

1* 




1 










D 

2» 


U 

1* 

1« 

1« 

u« 

c « 

c< 

1 « 



-- -i - - 

— 

■ 



- “ 


- 

- 



C 




c* 

1« 

0< 

1* 


1« 

1« 


. 

__ . - 

— 


■ - - -- 


— 

- • 

, - .... 









LI 

L2 

L3 

LA L 









L. 

n 



1 










e 

1* 

0* 

■ e*- 

- 0-»- 

1* 

2« 

■1^ 

2 ■» 

G« 

iM 


1 










c 

1 * 

0* 

o« 

0* 


2« 

1* 

2 « 

G « 

IM 


1 



















,n_*. 

— 4-^-- * 

B- 

r K 

— 


tj-H 

— t-R — 

JJ ^ 






1 










- -Q- 

. 1 - 


- fHt 

■ 0* 

1« 

2« 

1^* 

2 « 

CM 



1 










— — 'E- 

- - ■ 

-1. -» - 


— . 

. 


C* 

C « 

Qm 

iM 


1 










B- 

3-»_ 


— 6-»- 

— s-»— 

-Hh- 


— 1-»-. 

— 2-« 


l-M- 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 










Lr-t- * 

- L't 

tr 

• tr 1 

1? 


1 , 








2 

C 

i 


D 

i 


2 * 


c« 

1^* 

2» 

c» 

C « 

I « 

Hf 



c. 












D 


2« 


c* 


1* 

U 

1* 

C « 

1* 











LI 

L2 

L2 


L 




- 

- 

- 



3 

2 

1 

C 



1 












C 

2 

T «■ 

2h* 


0« 

u 

3« 

G-tt 

C « 

0 « 



0 


C» 

2* 


c* 

Hf 

3« 

C« 

C « 

■c« 

1* 



2 












- ••'C 

1 

- 

-1 If 


-2'*' 


1« ' 

1» 

i« 

■ e « 

1^* 


c 


3 * 

L * 


c* 


T^Tf — 

W ** 

— 2tt- 

— 1-«- 

— ct«- 

in- 



- - 




“ 




LL- ■ 

L2 

L3 


L 


_ 







■a 

2 

C 

c 


c 


C « 

3 « 

c« 

c« 

1« 

3« 

C« 

C « 

r, « 



D 

1 

3 » 

C* 

c* 

0* 

1« 

0« 

3 * 

C « 





L6 


16 


L6 



1 



-- 

- . , 

- - - 

- 


LI 

4 

L2 

Vp 

L3 

0 

L4 
2 - 

C 


2 « 


0 ■» 

I* ■ 


-2tt- 

-CTf- 

C'^'' 

i-Tt- 

-- iTT- 


1 











C 


1 « 

C « 


2« 

m 

3 « 



C^* 

i« 


1 











c 


2^* 

c»- 


1« 

i« 

2« 


c « 

1 « 

■ !« 


1 





iTf- 


-Ctt- 

-c-« 

- 2-^f 

'• H « 

L 


-?T5 

“tr*~ 

■■" C-» 

— 

— 


— 

- — 

— 

■ . 



L1 

L2 

L3 

L4 


“-1 








c 

1 

1 

C 

1 


c « 



i « 

3« 

C« 

c « 

C « 


C 


3 ■» 

c« 

c« 

c« 


1* 


1* 

1 
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Table 2-3. SERIES, FOR ORTHOGONALITY CONDITION (continued) 








■ 


LI 

L2 L3 LA 

Li 

L6 


1 







— ^ — 

— 1-' -C T 



c 

1 

1 * 

1 -B 

Q « 

1-B' 

1 « 

■3.5- 

— 

1 « 




2 

"TTf 

1 B 

0 * 

tTf 

-Itt- 



G-«- G'B” i'B 



u 

1 

-2ti- 

1 

C'B'“ 

0 * 

-tB- 

2«— 


-C“B— t-B- ■ —1 *- 



D 



C'»" 

0 b 

0 ^ 

-fir- 

— ty~ 

J-1T- 

— C"« — ir« tB- 

-L'i- 

•" t 6 









“l-t 

•L'Sf • 

11 C 

C 

T" 

1* 

1 « 

IB 

Cb 

IB 

3« 

0* 

C B C « 1 b 



D 


1» 

1 B 

IB 

OB 

1« 

3 B 

C-tt 

C'b c b 1b 



D 

1 ‘ 

2 » 

Cb 

Ob 

Cb 

IB 

IB 

L* 

1 B 0 * 1 « 











LI 

L2 L3 L^ 

Li 

L6 


1 







— A — 

— 2 g g 

' — -‘-1 •(( 



c 


■]7ir- 

2"«- 

— 0^ 

•■■'Cb - 


-5* - 













XT UZ 'trr- * L*1 

5 C C 1 

Li 

Lc 

c 


2* 

C * 

Gb 

IB 

IB 

3« 

Ob 

c B c * 1 B 



D 

"1- 

2 « 

C B 

Ob 

iB 

IB 

3 « 

Cb 

Ob C « 1 « 



C 

3 « 

c * 

Cb 

Ob 

iB 

2b 

Cb 

C b 1 b 1b 



■ 

2 

- 

- 


■ -G-B- 

iB 

3« 

L L 
5 

L2 L3 L^ 

i C 0 

• C'B c B- 1b 

Li 

L6 

c 


-2‘TT- 

• fis 

--■tj-v ■ 

0 « 

ti- 

- 1 1 - 


t- 







iri — 
6 

“It? “t“2“ ‘tr^ 

C C O' 

0 

1 

Cb 

c « 

G« 

IB 

3« 

• Cb 

C B C « 1 « 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 


Li L2 L3 L^ LS Lt 

; ^ ----- -• 

1 

■ C ' -O'* 2'« "0 1« im- - 2-^* C«'2« 

- • - Li L2 L3 LA Lf Lfc 

1 A I 1 

1 : 

C C« 2» C« I* 1* 2» 1 ■» 1« I* 

LI L2 L3 LA L5 L6 

-1 -A 2 G 

1 

— B C'-n 2 ^ Q-« 6^« 1-'« 2^* 2“« €•■«“ T-*- 

- - ----- ■ tl L2 L3 ■ LA LI Lt 

2 3 C 2 

C C* 2» 0« 0* 1« 2« 1« C« 2* 1* 


LL L2 L3 LA L3 Lfc 

- 2- -2-1 1 

1 

... — g. ^ 2-« — X -* — ■ 2r»- ■ i~* ■ -i* I* !•* 

t-i L-r-Hr2 IHr— - -L-fc- - - 

2 2 2 C 

1 - 

C C* 2« C^< 0« 1« 2» 1» 2« C« 

LI L2 L2 L* LS L6 

j Q -—2 — - -- 

1 • 

-- D " 2 ~*- e-^f- -<3-* --6-* - 1* -l•^t 2-»“ C« 2 * l« 

■ LL L2 L3 LA L5 L6 

3 2 11 

C 2* Q« C'« C« l-« 1* 2* i« 1« 1« 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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1« 

2» 

1* 
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L2 L3 

LA 
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L6 
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1. 
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LA 

L^ 
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C 
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2« 

1« 
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LA 

L5 
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C 
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i« 
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c 
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0« 
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1 « ■ 1 « 
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C 
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1« 
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1« 

1« 

2» 


1« 
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1 

C 
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c« 


1« 

2« 
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2« 


r 

c 

1« 
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0* 
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0« 

1« 
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c 

1« 

2« 

1* 

1« 

1« 

0« 

U 

1« 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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la 
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0 
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La 
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Table 2-3, SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES' FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3, SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGOMALITY CONDITION (continued) 
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Tabl© 2-3. SERIES FOR ORTHOGONALITY CONDITION (continuGd} 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (continued) 
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Table 2-3. SERIES FOR ORTHOGONALITY CONDITION (concluded) 
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SECTION III 

OPTIMAL RENDEZVOUS WITH SATELLITE IN CIRCULAR ORBIT 

3.1 INTRODUCTION 

The mathematical model for this problem is nearly the same as that for 
optimal ascent to circular orbit. The only difference is that one of the terminal 
conditions is different. Instead of a non- trivial transversality condition^ 
there is a condition that describes the rendezvous of the powered vehicle and 
a target satellite. 

The coordinate system is shown in the following figure. 



0 as the angular rate of the target 

(j, as the angle between launch vertical and the target at the time of 
"ignition". (Assumed to be known.) 
t as the time when the flight begins, the "ignition" time 

ig 

t as the current initial time 
o 

t^ as the time from t^ till rendezvous 
At as (tj - t^). 
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At any we compute (j) as 

♦ = *o + « <‘o - '^lg> ■ 

For the "rendezvous" terminal condition we have 


X. - R Sin + 0At) = 0 , 

f CO o ’ 

or - R Cos C A + 6At) = 0 . 

f CO 

In addition we have the terminal conditions for the circular orbit 


2 2 

YZ - V =0 

f CO 

2 2 „ 

R^ - R =0 

f CO 


R. 


Vj = 0. 


The solution of this problem proceeds in the same way as that of the optimal 
ascent to circular orbit. The solution for At that was obtained in that prob- 
lem applies here. 


The expression of the rendezvous condition is written as 


0 - ¥ + W,At + W^At + 

o 1 Z 


where 


¥ = X _ 


R sin + 6At) 
CO o 


¥, = u 


¥„ = u. 


¥., - 


u 


This expansion can be put in terms of the multipliers by substitution for the 
u time derivatives. Together with similar expansions for the final radius 
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and path angle conditions and the scaling conditions, we obtain a system of four 
algebraic equations in the multipliers. 


The unknown parameter, time-to-cutof f , or At, appears implicitly in the 
rendezvous condition. It was necessary to approximate the sine term with a 
series expansion in order to have At appear explicitly. The term is replaced 
by: 


W » X - R (d) + 9At - 0At 

o o co ^o o ^o 

- c|)^ 0^ At^ - 0^ At^). 


3.2 SOLUTION FOR THE MULTIPLIERS 


The solution for the multipliers in this problem was carried out in the 
same way as in the ascent-to-circular-orbit problem with the exception 
that the transversality I condition of the circular problem is replaced 
by the rendezvous condition. To obtain a reasonably accurate solution for the 
multipliers it was- necessary to have a fifth-order series representation for 
the rendezvous condition. The other conditions were represented by third-, 
fourth-, and fifth-order series. 

Typical solutions are shown in Table 3-1. The equations were solved xvith 
nominal values of At assumed. The Newton- Raphson method was used to compute 
solutions; three or four iterations were usually required before corrections 
became less than 0.5 x lO”'^. Initial estimates were the nominal multipliers' 
values . 

No reason was fotmd for the case where convergence did not occur. The 

-3 

tolerance on the corrections to successive Iterates was increased to 10 but 
there was no convergence. The derivations and coding for this system of equations 
is being checked for errors. 
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Table 3-1. SOLUTIONS OF EQUATIONS FOR MULTIPLIERS 


Nominal Multipliers 
0.97371974 
^2 = 0.22774957 
A 3 =^0.17809198 X 10'^ 
A^ = 0.45610352 x 10“^ 


Computed Solutions 
\^= 0.98382817 
A 2 = 0.17911486 
A 3 = -0.43654483 x 10 
A^ = 9.42815572 x 10“' 


Radius condition represented by third-order series in At 


Orthogonality condition represented by third— order series in At 
Rendezvous condition represented by fifth-order series in At 
At — 170.34 sec .Nominal Trajectory AA-1 


Nominal Multipliers 
A^ = 0.97371974 
X 2 = 0.22774957 
Xg =-0.17809198 x lO"^ 
= 0.45610352 x 10"^ 


Computed Solutions 
A^ = 0.96952223 
A 2 = 0.24500336 
A 3 =-0-. 17845442 x 10 
A^ = 0.48520438 x 10 


Radius condition represented by fourth-order series in a^ 


Orthogonality condition represented by third-order series in At 
Rendezvous condition represented by fifth-order series in At 
At = 170.34 sec Nominal Trajectory AA-1 


Nominal Multipliers 


Computed Solutions 
Non- Convergence 
Non-Convergence 
Non- Conver genc.e 
Non-Convergence 
Radius condition represented by fifth-order series in a t 
Orthogonality condition represented by fourth-border series in At • 
Rendezvous condition represented by fifth-order series in At 


Aj^ = 0.97371974 
A 2 = 0.22774957 
A^ =^0.17809198 x 10 
A^ = 0.45610352 x 10 


- 2 . 


-2 


At = 170.34 sec Nominal Trajectory AA‘-1 
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:>ECTI0N IV 

OPTIMAL ASCENT TO ELLIPTICAL ORBIT 

4.1 INTRODUCTION 

A third guidance problem studied was that of optimal ascent to elliptical 
orbit. The approach to solution is the same as for the previously described 
problems except that the terminal boundary conditions are different. It is 
assumed that the orbit is described by its associated, specific energy and 
magnitude of angular momentum, and the direction of perigee with respect to 
an earth- centered coordinate system. These terminal conditions may be 
expressed through the following equations: 

(Energy) 

(Momentiim) 

(Perigee) 

(Transversality) . 


^ ^ ^CO - » 


X 0 




[hi 


X. • + X. ’ = 0 

r r r r 


The first two conditions determine the major axis and eccentricity of the 
orbit. The third condition is derived from "Hamilton's Integral" for the two- 
body problem; <j> is the angle between the x-axis and the perigee direction. 

The coordinate system is shown in the following sketch. 


4-1 



NORTRONICS - HUNTSVILLE 


TR-792-8-283 



As noted, the approach to solution for this problem is the same as for the 
two problems described previously. A nominal trajectory (i.e., a numerical 
solution for optimal ascent to elliptical orbit) was obtained. Data from this 
nominal case was then used to test, numerically, the analytical developments. 
Attempts were made to achieve correlation between the nominal case and values 
predicted by analytical approximations. At the time that these comparisons 
were made, two-dimensional nominal trajectories were not available. A nominal 
case was used that was very slightly three-dimensional; out-of-plane flight was 
only a few meters. 

4.2 SOLUTION FOR FLIGHT TIME (At) 

In the optimal ascent to circular orbit and rendezvous problems, the ex- 
pression for At was obtained by inversion of one of the series expansxons for 
a terminal condition. This expression could then be substituted into the other 
expansions and the explicit appearance of At would be eliminated. However, it 
may be desired to solve five equations in five unknowns if At is not eliminated 

If At is eliminated by substitution, one should choose the "best" series 
to invert for At. The "best" series is the one that requires the fewest terms 
to implicitly define At and also requires the fewest terms to invert. To 
determine which series would be the best candidate for inversion, all four 
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series expansions were taken to third order and their coefficients were evaluated 
with nominal data. These third-degree polynomials in At were then solved for 
all of their roots. The third -order expansion for the terminal energy condition 
appears to define the most accurate values for At. 

In addition, approximations for At and its powers, analogous to those in 
Section II, were attempted. Let the series expansion be 


where 



etc. 


Asstime that an estimate to At, defined as x, is available. After following 
through the steps outlined in Section II, the following approximations are 
obtained: 


where 


At X + dt 


A 


2 2 

At = x'" + 2x dt. 


3 3 2 

At ~ X + 3x dt, , 
A 




W. ( 


T + A^T^ + A3T3) 


1 + 
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Figure 4-1 shows the result of using various values of t to estimate At and its 
powers. While the At values predicted by the formula above are not close to 
the true At values, they are close to the actual root of the polynomial that is 
closest to the true At value. If the series expansions of the terminal energy 
condition were taken to the fourth or fifth order, and then inverted, the 
approximations would probably be closer to the true values. 

However, because the approximations for At are rather complicated in this 
problem, it is advisable to consider solving for At simultaneously with the 
multipliers. There was not enough time available to follow through on this 
possibility during the term of this contract. 

4.3 SOLUTION FOR MULTIPLIERS 

In solving for the multipliers, it was asstimed that At would be eliminated 
by substitution of the inverted series for the energy condition. Thus, the 
series for the perigee, transversal ity, and' momentvun, together with the scaling 
condition on the multipliers, would be used to derive a system of algebraic equations 
that could be solved for the initial values of the multipliers. These algebraic 
equations were derived by substituting the differential equations' derivatives 
into the third-order expansions. 

Attempts to solve, numerically, this system of four equations in four un- 
knowns failed. Derivations and program coding were checked for errors but none 
were apparent. The expansion of the transversality condition to third order 
has a very large remainder term when evaluated with nominal data. It is likely 
that if this series is extended to higher orders, a solution for the multipliers 
can be achieved. Again, there was insufficient time available to pursue this 
problem during the tern of this contract. 
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Figure 4-1 shora the result ef using verlnus valnen of t to estiwte it and Its 
powers. Hhiln values predlotaa a» f«r«la above are nof close to 

the true 4t values, they are close to the actual root of the polyn«ial that la 
closest to the -tme « value. If the series expansions of the temlnal energy 
condltlnn nare tahsn to the fourth or fifth ordar, a»i then Inverted, the 
approximations mvU probably be closer to the true values. 

However, because the app»xl»tlons for 4t are rather cospHcated In this 
probles., it is advisable to consider solving for it siianltaneously with the 
BUltipllers. there was not enough time available to follow through on this 

possibility during th& term of tMs contract. 

4.3 SOLUTION FOR M0LTIFLIEES1 


In solving for the Mltlpllers, It was assumed that 4t would be eliminated 
by substitution of the inverted aeries for the energy condition. Ihns, the 
senes for the perigee, transversallty. and momentum, together with the sealing 
condition on tie mnltlpHer3,«onW he used to derive a system of algebraic equations 
that could be solved for the initial values of the multipliers, these algebraic 

equations were derived by substituting the differential equations' derivatives 

into the third-order expansions. 


Attempts to solve, numerloally, this system of four equations in four nn- 
Knoims failed. Derivations and program coding were checked for errors but none 
were apparent. lie expansion of the transversallty condition to third order 
has a very large remainder tea when evnlnated with nominal data. It is likely 
that if this series is extended to higher orders, a solutl-on for the multipliers 
can be achieved. Again, there ws Insufficient time available to pursue this 
probl«5 during tfee tarm of this contract. 
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SECTION V 

SYMBOLIC MANIPULATION OF POLYNOMIALS 


5.1 BACKGROUND 

An outstanding characteristic of the work described in this report is the 
amount of algebraic manipulations required. These manipulations can be put into 
two classes; (1) Differentiation, and (2) Substitution, multiplication, and 
simplification of polynomial-like expressions. 

There are presently available a mmiber of computer languages that are 
more or less suitable for handling these kinds of symbolic mathematical problems. 
One language that appears to be strongly oriented toward function differentiation 
is FORMAC. It, like the other symbol-manipulation languages known to the 
authors, suffers from an inability to exploit off-line on "non-core" storage. 

By the time the significantly difficult parts of our analytical develop- 
ments were reached, the storage facilities of FORMAC had been exhausted. 

Attempts were made to segment problems into portions that could be handled by 
FORMAC. But it was soon evident that the bookkeeping problems associated with 
this approach offset the advantages. This is not to say that FORMAC or a 
similar language is without value, but that for the problems we dealt with it 
was of relatively small utility. 

It was ap})arent that most of our symbolic computations involved polynomial- 
like structures; e.g., power series or determinant expansions. Most of the 
operations required with these expressions Involved multiplication as the most 
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difficult step; e.g., substitution of polynomials. These are also the kinds 
of operations that the symbol-manipulating languages fail to perform for 
polynomials that have a large number of terms. 

Because of these difficulties, it was decided to develop a computer 
program, written in a general purpose language, that would accomplish at least 
the multiplication of polynomials. 


5,2 THE MDLPO PROGRAM 

The MULPO Program is designed to multiply polynomials, in one to six 
variables, which have literal coefficients. An additional feature is the 
ability to add polynomials of the same type as it can multiply. Thus, the program 
can perform substitution of polynomials into other polynomials. It is written 
entirely in FORTRAN IV, and is intended for use on the IBM 7094. However, it 
can be run, with minor changes, on another machine with a comparable FORTRAN IV 
translator and enough tape units. A users guide for MULPO is given in reference 
3. 


At one point in its development, MOLPO was written in ALGOL 60 for 
the Burroughs B5500 computer. The reason for this was to obtain faster 
execution speed. Due to the widespread use of FORTRAN, however, it was con- 
verted to that language. For all but unusually long problems or volumes of 
production work, the relative inefficiency should not be noticeable. 

As a test of its ability to handle large-scale problems, the program was 
asked to multiply four polynomials together; then add this result to the 
product of three polynomials; then add this result to the product of two 
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other polynomials. Each polynomial was in terms of four variables and consisted 
of about ten terms each. The result was a polynomial of over 11,000 terms; and 
its like powers and coefficients were combined. 

In addition to studies such as reported here, a program like this 
probably has numerous other applications; for example, in perturbation techniques 
in celestial mechanics. 


5.3 DIRECTIONS FOR READING OUTPUT FROM MULPO 


Consider the following problem: 
Given : 


^o ^ -^ijkl ^1 
h = ^ «ijki 
h = ^ '=iaki 
h = ^ "ijki 
^ \jki 

J \jkl 


x-^ X^ X^ 
^2 ^3 ^4 


2 3 4 


xj X^ X^ 

^2 ^3 ^ 


xJ 

^2 ^3 ^4 


x3 X^ 
^2 ^3 \ 


2 3 4 


^3 = I Z., 


x"^ X^ x^ x^ 
ijkl ^1 ^2 ^3 ^^4 


where 


\jkl ^ ^ijkl 

"ijkl = ^ ^ijkl 


Form 


S = ^ S 


ijkl '^l 




^1 = 


■’o + 


+ P3Q, 
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and Simplify” s by collecting lilce combinations of powers resulting from 
multiplications, and by combining like coefficients when they occur as sums 
in a like combination of powers. For this example, aii S-polynomial term will 
generally have a multi-term coefficient made up of a term from X^-eoefflcient)* 
a term from the product of (B-coefflcients times X-coeff icients) , a term 

from ^2^2 Icients times Y-coefficients) , and a term from (D- ' 

coefficients times 2-coefficients) . - 


Referring to Table 2-2, it can be seen that the integers directly under 
LI, L2, etc are the powers of the lambdas in the S~polynomial. (Note that the 
output is arranged so that the polynomial is ordered according to increasing 
powers.) The integer on the next line is the numerical multiplier of the 
literal coefficient oh the next line down. 

The symbol for the coefficient comes first, then the next four integers 
represent its subscript. The fifth integer is the exponent on the coefficient. 
(The asterisks should be ignored when they, appear . ) When the remainder of the 
line is blank, interpret this as a plus sign and go to the next line. If- only 
one more term appears (i.e., five integers), then these Integers represent the 
four subscripts and the exponent of the coefficient in a Q-poiynomial that 
multiplies the term immediately to the left. 


Which Q-pol 3 momial is involved depends on the coefficient symbol at the 
start of the line. The coefficients of the polynomial- always multiply the 
B-coefficient, the polynomial's coefficient the C-coeff icients , etc. 


5-4 



NDRTRONICS - HUNTSVILLE 


TR-792-8-283 


If more than three terms are on one line, the last term Is a coefficient 
of a Q-polynomial. The preceding coefficients are, as before, identified by the 
alphabetic symbol at the far left. In some cases enough coefficients are 
combined to cause the program to carry over the remainder of a coefficient to the 
next line. In such cases the appropriate alphabetic symbol reappears at the 
beginning of the next line. 


Actually, reading the output is far simpler than describing how to do so. 
We can translate the first four terms of Table 2-2 into conventional notation 
as : 


A + 
^0000 

®0000 ^0000 
^0000 "^0000 
®0000 ^0000 

^0001 ^0000 

^0010 ^0000 

®0000 ^0100 
Sooo ^0100 
^0100 \iooo 
®0000 ^0100 
°0100 ^0000 


Coefficient of X 


0 

1 


x° x° 


X 


0 

4 


Coefficient of X^ X^ X^ X^ 


Coefficient of X^ X^ X^ X^ 


Coefficient of X^ X^ X^ X^ 
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SECTION VI 

CONCLUSIONS AND RECOMMENDATIONS 

Based on tlie studies conducted under this contract, several conclusions 
and recoiBateadations can be made. The first part of this section deals with the 
study described in this report. The second part is concerned with recommended 
modifications to the approach described herein and possible alternative 
approaches to obtain explicit guidance formulas. 

A major accomplishment has been the successful automation, on digital 
computers, of the manipulation of polynomials that have literal coefficients. 

This computer program has an advantage over other, general purpose symbol 
manipulating languages because there are no practical limitations on computer 
storage. In conjunction with a langttage suchas foR^G, it should give an 
analyst the capability of multiplying or substituting power series or poly- 
nomials of virtually any number of terms. The program was designed to accommodate 
polynomxals or power series involving as many as six variables, should the need 
arise in the future. It was also coded in a general purpose language (FORTRAN 
IV) that is commonly used and should be available in the future. It is also 
possible, but not easy, to alter the program to perform partial differentiation 
of multivariable power series or polynomials. 

Much of the effort during this study was directed toward simplification 
of the analytical developments. The term "simplification'* is taken to mean 
the dropping of terms in formulas for the guidance functions when these terms 
do not significantly influence the answers, or approximating complicated 
formulas with simpler ones. The simplification process was rather successful 
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in the formulas for At; particularly for ascent to circular orbit and rendezvous 
in circular orbit. However, no simplifications in the solution for the 
multipliers or guidance function beyond those previously reported were achieved. 
A further test of these simplifications was made by perturbing the initial 
values of the state variables on a nominal trajectory. Simplifications that 
were valid on the nominal trajectory were also valid on neighboring trajectories. 

Application of the Successive Substitution method to solve the equations 
in the multipliers is recommended only when the coefficients of the system of 
equations are numbers. Application to a system having literal coefficients 
results in an impractical number of terms. If approximate values of the initial 
multipliers are assumed known, the method of inversion should be tried. This 
approach to solution of the equations in the multipliers is analogous to that 
used in Section II for solution for At. Such an approach could allow for a 
considerable amount of simplification as was found in the approximations for 
At. 


No solution was obtained for the multipliers in the problem of optimal 
ascent to elliptical orbit. This is believed to be due to the insufficient 
order of the series expansion for the transversality condition. It is 
recommended that this problem be studied further, series orders extended, and 
solutions for the multipliers be attempted. Because of the complexity of the 
formula for At that is got by inversion of the energy condition expansion, 
consideration should be given to solving for At simultaneously with the 
multipliers. 
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SECTION VII 

BOLZA PROBLEMS WITH END ORBITS 

The work discussed in reference 2 pertaining to Bolza problems is not 
closely related to the topics discussed in other parts of this Summary Report. 
However to keep the reader informed on all aspects of the contract efforts, an 
abstract of reference 2 is provided in the following paragraphs. 

It is well known that many of the problems of space vehicle guidance can 
be formulated in such a way as to be susceptible to the methods of the calculus 
of variations. The problems considered fall roughly within the class of problems 
referred to as problems of Bolza. 

Many space missions are formulated in such a way that considerable further 
analysis is needed to bring the problem within the format of the Bolza problem, 
as regards the end-point constraints. It is usually the case that conditions 
are imposed on the orbits determined by the end points of the burn trajectory, 
but the translation of these conditions into conditions applying directly to the 
end point values of the states and time requires that the orbital motion be 
expressible as a function of the end-point conditions and time, and that the 
conditions on that orbital motion then be invertible to yield conditions on the- 
end-point values. Both of these problems can be quite profound, so that it is 
not often possible to give analytic expressions for the end-point constraints , 
although „ these are assumed available in the formulation of the problem of 
Bolza. 
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We therefore formulate a modified Bolaa problem, called a Bolza .problem 
:«ith_enl_or^. The objective is to define the problem in terms of that 
xnformation which actually is available, and to obtain conditions which apply 
to that available information. Thus in the new problem formulation, there 
are given differential equations (analogous to the gravity equations) to which 
the end points of any admissible curve . determine particular solutions, called 
the end orbits of that curve. Other conditions are given, which apply to the 
end orbits (uot to the end points), called end-orbit constraints, or simply 
orbital constraints. The transversality condition is then deduced in such a 
form as to apply directly to the differential equations of the end orbits and 
the orbital constraints. In particular, it is not required that the differential 
equations of the end orbits be solved. 
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APPENDIX A 

TIME DERIVATIVES OF DIFFERENTIAL EQUATIONS 

V = I - 

— * * 

V = a^A + a^x - a^(I • X)I - - 3^V 

V = a^X + “ 3a2(X • X)X + A^(T • R) 

• « • • • 

-2a^(X • X)X - a^cT • xjT 4- a^(X • T)^ T + 

X = g^X + y(X • R)R 

X = e^X + 8^X + |^p(X R) + y(X • R + X • X)j R +y(X * R)V 

where 


a = £ o, =_£-S „ _ Z H 

1 Ki’ 2 in m* 3 m m 


a , GM 

“ 3 » ~5 

R ^ \r^/ 


(R • V) 


A, = Ca,Y + a^6^)R, y = — ,P = 15 |p| (R • V) 


^2 " [4 ■” * V) + g^(V • V)] R - 262 V 


All variables are evaluated at t « t . 

o 
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APPENDIX B 

INVERSION OF SERIES IN SEVERAL VARIABLES 

Suppose that we are given the two equations in two variables. 


’*2^ = 0 


F2<Xj,*2) = 0. 


We wish to solve them for roots and x^. Let the functions and F 2 
be expanded about a point x^ and X 2 which is presumably near Xj^ and X 2 * 

Fi(Xi,X2) - ~ ^^2 ' * 2 ^ 


= 0 


and similarly for F 2 * These expansions are rewritten as 


“00 “10 *1 “01 ’'z “20 ’ 


600 = SlO + »01 =^2 »20 ='1 + "■ 


where 


^1 === ^^1 - ^1^ 

X2 = (^2 - X2) 

-«00 = ^1<^1’^2^ 
-8qo = F2^^1’^2^* 


We assume that the inverse ma.y be represented as 

^1 ^ ^10 “00 ^01 ^00 ^20 “00 ^11 “00 ^00 ^02 ^00 ’’’ 

^2 ^ ®10 “00 ®01 ^00 ®20 “00 \l “00 ^00 ®02 ^00 


B-1 



NORTRONICS - HUNTSVILLE 


TR>?92-8~283 


It is tiow necessary to deterniine the coefficients of the Inverse series in 
terra® of the coefficients in the original series. The derivations tfill be 
simpler if -we write the inverse as 


+ Z2 Zg + « ♦ • 

+W2+W3 + -, 

where Z^ and stand for all terms of degree n in the a^s and gfs. 


The assumed form for the inverse is substituted into the original series 

and the method of undetermined coefficients Is used to determine the Z and 

n 

first-degree terms we have 




“10 

%1 

-1 

**00 

^1 


®10 

^01 


^00 


For the second-degree terms we have 


Z2 


®10 

«01 

-1 


^2 


^10 

^01 


^1 


■% “20 ^ “11 h “02 ®i 

-Jj = 630 zl + Sn Zj Wj + S02 wj 


In general, 
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APPENDIX C 

SIMPLIFIED SERIES FOR At FOR ASCENT TO CIRCULAR 
ORBIT AND RENDEZVOUS PROBLEMS 


At = 


where 


0000 


1000 


^0000 ■*' %000 ^1 ^ ^0100 ^2 101.^1 ^2 H ®0201 4 ^4 

^0120 ^2 ^3 + ®1002 ^1 ^4 + ^SOOp %00 ^2 

^2111 ^2 S ^4 ®1202 4 4 ^ ®2120 ^1 ^2 ^3 

■ ^1211 ^1 ^2 ^3 ^4 *'^ ^1011 ^3 ^4 ^0111 ^2 ^3 ^4 

^220 S 4 4 ^ ^2102 4 h 4 ^ ^0210 4 h ^ ^2001 4 H 


[6lV^ + 32^^ * V)- aj - r^/2S + t + Z/(l - 2C) 

-[302V^ - * V) + p(R . v)2 + 

" UgU - 2ajBj^x] r^/28 


[a^u - 3a2 02^x - 3p20!j^x - y^u - 2a yxyv] x^/66 


^0100 ~*-“2^ ^ ^ T^/25 

la^v - y ~ SB^ctj^y - 4a^B^v - a^yx^v - Za^yxyu] r^/66 

^1101 ~ ~ [2a2U - t^/65 

^0201 ~ ~ f«3^v]T^/25- t2«2V - 3aj-B^y1 t^/66 


^0120 = ^“1^3x3/66 
^1002 ~ /6B 

^3000 ^ " ctj^Yuy^ - 2a^Yxyv]T^/65 
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®0300 ” ■ 20jTi!iyu]T^/66 

®21U - - 

®1202 ^ ■ I30ju]t^/6S 
*2120 = - 

*1211 “ ■ f7ajvlT*/6« 

*1011 - [2V’'"''** 

*0111 “ [2Vl'"/«* 

*1220 * PVI'"/** 

^2102 ^ t3a^v]T^/6S 

^0210 ~ / 2 ^ + 2a^u - 3 r^/66 

®2001 ~ ^/26 + 2a^v - 3 B^a^y] r^/66 

where 

== GM/R^ 

&2 = - 3(QJ/R^)(R-V) 

Y = -3GM/R^ 

P = 15 (CM/R^)(R‘V) 

6 ~ (v4) + (^4) 

Oj^ = F/m 

ct2 = -F/m (m/m) 

JL * * 

G = V . V/2V ' . V 

■*[ = the estimated value of At. 
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